Light scattering detection of quantum phases of ultracold atoms in optical lattices 
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Ultracold atoms loaded on optical lattices can provide unprecedented experimental systems for 
the quantum simulations and manipulations of many quantum phases. However, so far, how to 
detect these quantum phases effectively remains an outstanding challenge. Here, we show that 
the optical Bragg scattering of cold atoms loaded on optical lattices can be used to detect many 
quantum phases which include not only the conventional superfluid and Mott insulating phases, 
but also other important phases such as various kinds of density waves (CDW), valence bond solids 
(VBS), CDW supersolids and VBS supersolids. 
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Various kinds of strongly correlated quantum phases 
of matter may have wide applications in quantum infor- 
mation processing, storage and communications It 
was widely believed and also partially established that 
due to the tremendous tunability of all the parameters 
in this system, ultracold atoms loaded on optical lat- 
tices (OL) can provide an unprecedented experimental 
systems for the quantum simulations and manipulations 
of these quantum phases and quantum phase transitions 
between these phases. For example, Mott and superfluid 
phases may have been successfully simulated and ma- 
nipulated by ultra-cold atoms loaded in a cubic optical 
lattice |3j. However, there are still at least two outstand- 
ing problems remaining. The first is how to realize many 
important quantum phases [![. The second is that as- 
suming the favorable conditions to realize these quantum 
phases are indeed achieved in experiments, how to detect 
them without ambiguity. In this paper, we will focus on 
the second question. So far the experimental way to de- 
tect these quantum phases is mainly through the time of 
flight (TOF) measurement Q, H] which simply opens the 
trap and turn off the optical lattice and let the trapped 
atoms expand and interfere, then take the image. The 
atom Bragg spectroscopy is based on stimulated mat- 
ter waves scattering by two incident laser pulses 0, 5| 
through the TOF measurements. The momentum 4\ 
transfer Bragg spectroscopy was used to detect the Bo- 
goliubov mode inside an BEC condensate. The energy 
transfer Q Bragg spectroscopy was used to detect the 
Mott gap in a Mott state in an optical lattice. Opti- 
cal Bragg scattering (Fig.l) has been used previously to 
study periodic lattice structures of cold atoms loaded on 
optical lattices It was also proposed as an effective 
method for the thermometry of fermions in an optical 
lattice and to detect putative anti-ferromagnetic (AF) 
ground state of fermions in OL [U . There are very recent 
optical Bragg scattering experimental data from a Mott 
state, a BEC and AF state [9J. The atom Bragg spec- 
troscopy and Optical Bragg scattering are two different, 
but complementary experimental methods. 



In this paper, we will develop a systematic theory of 
using the optical Bragg scattering ( Fig.l) to detect the 
nature of quantum phases of interacting bosons loaded in 
optical lattices. We show that the optical Bragg scatter- 
ing not only couples to the density order parameter, but 
also the valence bond order parameter due to the hop- 
ping of the bosons on the lattice. At integer fillings, when 
q matches a reciprocal lattice vector K of the underlying 
OL, there is an increase in the optical scattering cross 
section as the system evolves from the Mott to the SF 
state due to the increase of hopping in the SF state. At 
1/2 filling, in the CDW state, when q matches the CDW 
ordering wavevector Q n and K , there is a diffraction peak 
proportional to the CDW order parameter squared and 
the density squared respectively (Fig. 3a), the ratio of the 
two peaks is a good measure of the CDW order param- 
eter. In the VBS state, when q matches the VBS or- 
dering wavevector Qk, there is a much smaller, but de- 
tectable diffraction peak proportional to the VBS order 
parameter squared, when it matches K, there is also a 
diffraction peak proportional to the uniform density in 
the VBS state (Fig. 3b). All the diffraction peaks scale as 
the square of the numbers of atoms inside the trap. All 
these characteristics can determine uniquely CDW and 
VBS state at 1/2 filling and the corresponding CDW su- 
persolid and VBS supersolid slightly away from the 1/2 
filling. In the following, we just take 2d optical lattices 
as examples. The Id and 3d cases can be similarly dis- 
cussed. 

The Extended Boson Hubbard Model (EBHM) with 
various kinds of interactions, on all kinds of lattices and 
at different filling factors is described by the following 
Hamiltonian [i,[lM3|: 

H B h = -ty^jblbj + h.c.) — fx ^ rij + — ^ m (m — 1) 
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where = b\bi is the boson density, t is the nearest 
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FIG. 1: Optical Bragg scattering of cold atoms moving in 2 
dimensional optical lattices. The q = ki — ki and uj = oji — ui2 
are momentum and energy transfer from the laser beams to 
the cold atoms respectively. The A stands for a aperture, the 
PM stands for a Photomultiplier. The off resonant scattering 
processes lead to the on-site term (a) and the off-site term (b) 
in Eqn.3. 
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FIG. 2: The charge density wave (CDW) phase in a square 
lattice at no = 1/2 with ordering wavevector Q n = (n, n). (b) 
valence bond solid (VBS) phases with ordering wavevector 
Qk = (tt,0) where the kinetic energy (Kij) = (b]bj + h.c.) 
takes a non-zero constant K in the two sites connected with 
a dimer, but in the two sites without a dimer. (c) Stripe 
CDW order at Q n = (tt, 0) and (d) Plaquette VBS order at 
Q n = (7r,0),(0,7r). 0113. 



neighbor hopping which can be tuned by the depth of the 
optical lattice potential, the U,Vi,V2 are onsite, nearest 
neighbor (nn) and next nearest neighbor (nnn) interac- 
tions respectively, the • • • may include further neighbor- 
interactions and possible ring-exchange interactions. The 
filling factor n — N a /N where N a is the number of atoms 
and ./V is the number of lattice sites. The on-site inter- 
action U can be tuned by the Feshbach resonance Q. 
Various kinds of optical lattices such as honeycomb, tri- 
angular [lH, body-centered-cubic 18 1, Kagome lattices 
[19( can be realized by suitably choosing the geometry 
of the laser beams forming the optical lattices. There 
are many possible ways to generate longer range interac- 
tion V\, V2, ■■■■ of ultra-cold atoms loaded in optical lat- 
tices. Being magnetically or electrically polarized, the 



or polar molecules |21| m K + 87 Rb ( 



52 Cr atoms [2 

or 39 K + 87 Rb ) interact with each other via long-rang 
anisotropic dipole-dipole interactions. Loading the 52 CY 
or the polar molecules on a 2d optical lattice with the 
dipole moments perpendicular to the trapping plane can 
be mapped to EqnQ] with long-range repulsive interac- 
tions ~ p 2 /r 3 where p is the dipole moment. The CDW 
supersolid phases studied by QMC [ll| and described in 
[15| by the dual vortex method was numerically found to 
be stable in large parameter regimes in this system 22 1. 
The generation of the ring exchange interaction has been 
discussed in 24(. Some of the important phases with 



long range interactions are listed in Fig. 2. Recently, the 
quantum ent ang lement properties of the VB state was 
addressed in [25| |. 

The interaction between the two laser beams in Fig.l 
with the two level bosonic atoms is: 

H int = [d 2 r*t(r)[ 1 f- + VoL(f f ) + ^a z 
J 2m a 2 



n 



(2) 



where *f/(r) = (ipc^g) is the two component boson an- 
nihilation operator, the incident and scattered lights in 
Fig. la and the two incident lights in Fig. lb have frequen- 



cies uji and mode functions ui{r) — e jfci ' r+i ^ ! . The Rabi 
frequencies VL are much weaker than the laser beams ( not 
shown in Fig.l ) which form the optical lattices. When 
it is far off the resonance, the laser light-atom detunings 
A; = oji — ui a where uj a is the two level energy differ- 
ence are much larger than the Rabi frequency and the 
energy transfer uj = cji — 0J2 ( See Fig. la and lb ), so 
Ai ~ A2 = A. After adiabatically eliminating the up- 
per level e of the two level atoms, expanding the ground 
state atom field operator if> g (r) = biw(f— fty in Eqnj2] 
where w(f — f\) is the localized Wannier functions of 
the lowest Bloch band corresponding to Vol (7) and 6, 
is the annihilation operator of an atom at the site i in 
the EqnlTJ then we get the effective interaction between 
the off-resonant laser beams and the ground level g: 
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where the interacting matrix clement is J^j = J drw(f— 

Ti)u\{r)u2{r)w(f — fj) = Jj t i. The first term in Eqn(3] 

is the on-site term D = Ji,i n i ( See Fig. la ). The 

second term is the off-site term ( See Fig. lb ). Because 

the Wannier wavefunction w(r) can be taken as real in 

the lowest Bloch band, the off-site term can be written 

as K = T,<ij> J i,M b j = J2<ij> J iA b \ h i + h - c -) which is 
nothing but the off-site coupling to the nearest neighbor 

kinetic energy of the bosons K%j = b\bj + h.c 

It is easy to show that: 



JV 



D(q) = fo(q)J2< 



' l m = Nf (q)n(q) (4) 



where q = k± — fcg, fo(o) — J dfe lq r w 2 {r) and n(q) 



Er bib? , - is the Fourier transform 
k f. k+q 



of the density operator at the momentum q. Note that 
n(q) = n(q + K). The wavevector is confined to L _1 < 
q < a^ 1 where the trap size L ~ 100/xm and the lattice 
constant a ~ 0.5/irn in Fig.l. In fact, more information 
is encoded in the off-site kinetic coupling in Eqn|3J In a 
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square lattice, the bonds are either oriented along the x matic increase of the scattering cross section 
axis fj — fi = x or along the y axis fj — fi — y, we have: 

K U = N[f x (q)K x (q) + f y (q)K y (q)} (5) 

where K a (q) = ± J^iLi er %9 ' ?i K i>i+a = 
e lq °</ 2 yZt cos kabtbr, - are the Fourier transform of the 

— ™ k K ~rQ 

kinetic energy operator K%j = b\bj + h.c. along a = x,y 
bonds at the momentum q and the "form" factors f a (q) = 
f(q,fi — fj = a) = J dfe~ lq ' r w(f)w(f + fi — fj). Note 
that K a (q) = K a (q + K). Following the harmonic ap- 
proximation used in [2j, we can estimate that /o(7r, 0) ~ 

e-Hyo/Br)- 1 " J x fa ) „ le -HVo/E r )- 1/2 -4(Vo/E r )V^ 

so |/x(7r,0)//o(7r,0)| ~ e ~T^V^ where V and 
E r = h 2 k 2 /2m are the strength of the optical lattice 
potential and the recoil energy respectively |2|. The 
/o(7r, 0) is close to 1 when Vo/E r > 4. It is instructive to 
relate this ratio to that of the hopping t over the onsite 
interaction U in the Eqn[TJ \f x (ir, 0)// (tt, 0)| ~ 
where a s is the zero field scattering length and 
a = A/2 = 7r/fc is the lattice constant, using the typical 



(7) 



across the Mott to the SF transition due to the pref- 
actor TV 2 . This prediction could be tested immediately. 
Surprisingly, there is no such optical Bragg scattering 
experiment in the superfluid yet. 

In the CDW with Q n — (ir, tt) in Fig. 2a, due to the lack 
of VBS order on both sides, the second term in EqnJS]can 
be neglected, so that 



da 
dVLdE 



\CDW 



& 2 iVU(<f)| 2 M<7>) (8) 



values t/U ~ 10 -1 , a s /a ~ 1CP^, one can estimate 
\fa/fo\ ~ 10 -3 . Note that the harmonic approximation 
works well only in a very deep optical lattice Vq 3> E r , 
so the above value underestimates the ratio, so we 
expect |/ a // | > 10- 3 . 

The differential scattering cross section of the light 
from the cold atom systems in the Fig.l can be calcu- 
lated by using the standard linear response theory: 



da 
dndE 



S(q,u) 



- (^) 2 ^ 2 [|/o(g1l 2 S„(g>) 
+ E \f«®\*SKM,<»)] (6) 



where q = k% — ko,cu = u)\ — ui2, the S n (q,u>) = 
(n(—q,~uj)n(q,uj)) is the dynamic density-density re- 
sponse function whose Lehmann representation was 
listed in 0]. The S Ka (q,u) = (K a (-q, -u)K a (q, u)) 
is the bond-bond response function whose Lehmann rep- 
resentation can be got from that of the S n (q,Lu) sim- 
ply by replacing the density operator n(q) by the bond 
operator K a (q). The integrated scattering cross sec- 
tion over the final energy = J dE (j fZ E is propor- 
tional to the equal-time response function 4£ = S(q) ~ 

(^) 2 N 2 [\f (q)\ 2 S n (q)+E a=£ , y \Uq)\ 2 SKA<l)}- 

We first look at the superfluid to Mott transition at 
integer filling factor n. When q is equal to the short- 
est reciprocal lattice vector K = (2%, 0), in the Mott 
state, ^ - |/ M (2vr,0)| 2 iV 2 n 2 , in the superfluid state, 

^f- ~ |/ ^(27r,0)| 2 7V 2 n 2 + 2|/^(27r,0)| 2 7V 2 J B 2 where 
B is the average kinetic energy on a bond in the super- 
fluid side. Because |/ sf (2tt, 0)| 2 - |/ m (2tt, 0)| 2 - 1 and 
B is appreciable in the superfluid side, we expect a dra- 



which should show a peak at q = Q n ( Fig. 3a ) whose 
amplitude scales as the square of the number of atoms 
inside the trap ~ |/o(7r, 7r)| 2 7V 2 m 2 where m — tia — ns 
is the CDW order parameter [15]. When q = K, 
then Scdw(K) - |/o(2tt, 0)| 2 iV 2 n 2 where / (27r,0) - 
/q(7t,7t) ( Fig. 3a ). So the ratio of the two peaks in 
Fig. 3a is ~ m 2 /n 2 if one neglects the very small difference 
of the two form factors. Slightly away from 1/2 filling, 
the CDW in Fig. 2a may turn into the CDW supersolid ( 
CDW-SS ) phase through a second order phase transition 
15J. Then we have (n(q)) — m5~^ + nd^Q where n = 
tia + rig = 1/2 + Sn. The superfluid density p s ~ Sn — 
n — 1/2. The scattering cross section inside the CDW- 
SS: ScDW~Ss{Qn) ~ liol 71 ", 7r)| 2 7V 2 m 2 stays more or less 
the same as that inside the CDW, but Scdw-ss(K) ~ 
\fo(2Tr,0)\ 2 N 2 n 2 + 2\f x {2Tr,0)\ 2 N 2 {6n) 2 B 2 will increase. 
The B is the average bond strength due to very small 
superfluid component p s ^ Sn — n~ 1/2 flowing through 
the whole lattice. So the right peak in Fig. 3a will in- 
crease due to the increase of the total density and the 
superfluid component inside the CDW-SS phase. 

Now we discuss the VBS state with Qk — (tt, 0) in 
Fig. 2b. Due to the uniform distribution of the density 
in the VBS, when q — K, the second term in EqnjS] 
can be neglected, so there is a diffraction peak ( Fig. 3b 
) whose amplitude scales as the square of the num- 
ber of atoms inside the trap ~ |/o(27r, 0)| 2 V 2 n 2 where 
fo{2n, 0) ~ /o(T)0) and n = 1/2 is the uniform density 
in the VBS state. However, when one tunes q near Qk, 
the first term in Eqn|6]can be neglected, then 

a—x^y 

which should show a peak at q — Qk signifying the 
VBS ordering at Qk whose amplitude scales also as 
the square of the number of atoms inside the trap ~ 
\f x {ir,0)\ 2 N 2 K 2 where K = K x — K y is the VBS or- 
der parameter [l5j]. So the ratio of the VBS peak at 
q = Qk over the uniform density peak at q — K is 
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FIG. 3: The optical scattering cross section in (a) CDW, 
the ratio of the peak at Q n over that at K is ~ m 2 /n 2 ~ 1. 
(b) VBS state, the ratio of the peak at Qk over that at K 
is ~ \fx/fo\ 2 K 2 /n 2 > 10" 5 , but still should be visible in the 
current optical Bragg scattering experiments. 



- K 2 /n 2 \f x (Tr,0)/f (2Tr,0)\ 2 > 1CT 5 . However, the 
smallness of \f x \ 2 is compensated by the large number 
of atoms N ~ 10 6 , \f x \ 2 N 2 = (\f x \ 2 N) x N ~ N ~ 10 6 . 
Therefore, the Bragg scattering cross section from the 
VBS order is > 10" 5 smaller than that at q = K at 
the same incident energy ij„ ( Fig. 3b ), but still ~ 10 6 
above the background, so very much visible in the cur- 
rent optical Bragg scattering experiments. Slightly away 
from 1/2 filling, the VBS may turn into VB Supersolid 
(VB-SS) through a second order transition (Isj . We have 
(K x (q)) = BS^o + K5tf K and (n(q)) = (Sn + 1/2)S^ . 
The superfluid density p s ~ Sn = n — 1/2. The 
scattering cross section inside VB-SS: Svb-ss(Qk) ~ 
\f x (ir, 0)\ 2 N 2 K 2 stays more or less the same as that in- 
side the VBS, but S VB -ss(K) ~ |/ (2tt, 0)| 2 Y 2 n 2 + 
\f x (27T,0)\ 2 N 2 (Sn) 2 B 2 x + \f y (27r,0)\ 2 N 2 (Sn) 2 B 2 where 
n = 1/2 + Sn and the B x ,B y are the average bond 
strengths along x and y due to very small superfluid com- 
ponent p s ~ Sn — n — 1/2 flowing through the whole 
lattice. So the right peak in Fig. 3b will increase due to 
the increase of the total density and the superfluid com- 
ponent inside the VB-SS phase. Very similarly, one can 
discuss the VBS order at q = Qk = (0,7r). For the pla- 
quette VBS order in Fig. 2d, then one should be able to 
see the Sk(q) peaks at both (n,0) and (0,7r). So the 
dimer VBS and the plaquette VBS can also be distin- 
guished by the optical Bragg scattering. 

In this paper, we only focused on the optical Bragg 
scattering detections of the various ground states in a 
square lattice. The detections of the excitation spectra, 
the generalization to frustrated lattices, the effects of fi- 
nite temperature and a harmonic trap will be discussed 
in a future publication. 
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